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VEHICLE ATTITUDE DETERMINATION WITH 
A SINGLE ONBOARD CAMERA 

By Richard J. Bendura, Allen B. Henning, 
and Robert E. Smith, Jr. 

Langley Research Center 

SUMMARY 

Application of a postflight photographic method for determining the attitude history 
of a vehicle equipped with a single onboard camera is presented. Also included is an 
appendix showing the equation and geometric relationships of the method. For applica- 
tion, a minimum of two earth- related points need be identified on each photograph and 
the location of the vehicle with respect to the earth must be determined. In addition, 
the orientation of the camera axes system with respect to the vehicle is required. The 
method requires no other onboard instrumentation. 

The method was applied in determining the spacecraft attitude for a portion of the 
first flight test of the Planetary Entry Parachute Program (PEPP). The spacecraft was 
equipped with two relatively lightweight and inexpensive 16-mm motion-picture cameras 
which photographed the surface of the earth. Data were obtained for altitudes ranging 
from 130 000 to 143 000 feet (40 to 44 kilometers). A comparison is made of the results 
independently obtained from each camera with ground-based tracking data and with an 
azimuth datum point from an onboard magnetometer direction finding system. 

INTRODUCTION 

For many flight projects, data would be more meaningful if a time history of vehi- 
cle attitude could be obtained. Various instruments such as accelerometers, gyro plat- 
forms, rate gyros, horizon sensors, and cameras have been used singly and in combina- 
tion to define vehicle attitude. Of the cameras used, most have been restricted to highly 
sophisticated and expensive systems which are capable of obtaining very accurate 
results. However, for some flight applications such accuracy is not necessary. Where 
vehicle orientation angles within several degrees are sufficient, photographs from an 
inexpensive camera not designed as a metric instrument may be used. In reference 1 a 
method of vehicle attitude determination with a single onboard camera is shown for a 
special case in which the horizon and one identifiable landmark are available on each 



photograph. References 2 and 3 present an analysis requiring a minimum of four land- 
marks assumed to be at the same elevation in order to determine vehicle attitude. 

The data presented in this report were obtained by means of 16-mm motion-picture 
cameras which were intended only to photograph the deployment of an experimental para- 
chute (refs. 4 and 5). However, during analysis of the data, it was noticed that through- 
out a large portion of the flight several landmarks could be identified on the film. This 
provided an opportunity for determining the spacecraft attitude history - in this case 
through application of a camera orientation method. The appendix presents the equations 
for the camera orientation method used in this investigation and is intended to serve as a 
guide for future applications. The method requires a minimum of two landmarks within 
view of the camera, location of the camera with respect to the landmarks, and orienta- 
tion of the vehicle with respect to the camera. 

SYMBOLS 


c focal length of camera 

h height above sea level of camera, feet (meters) 

I point in image space 

K constant of proportionality between object space coordinate system and 

image space coordinate system 

n number of observations 

Q point in object space 

R geocentric earth radius, feet (meters) 

t flight time, seconds 

X,Y,Z object space coordinate axes with origin at focal point; also, distance along 
these axes, feet (meters) 

X,Y,Z distance from geocentric coordinate axes origin, feet (meters) 

X e ,Y e ,Z e earth-fixed axes system 
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XpYpZj 


x,y 




r y 

Aff 

e 


image space coordinate axes with origin at focal point 

measured image coordinates in image space coordinate system 

coordinates of point of intersection between focal plane and focal axes in 
image space coordinate system 

spacecraft flight-path angle measured in Xg,Zg plane, degrees 
spacecraft flight-path angle measured in x E» Y E plane, degrees 
change in attitude angle 
error in x 


error in y 


A 


longitude of camera measured positive east from Greenwich, degrees 


X 3 ,m 3 ,i / 3 


J th 


directional cosines of image space coordinate system relative to object 
space coordinate system 

attitude angles with random error for error analysis, degrees 
theoretical attitude angles for error analysis, degrees 


OjjOg, <*3 camera azimuth, pitch, and roll Euler angles relative to earth-fixed axes, 

degrees 


<f> geodetic latitude of camera measured positive north, degrees 

$>' geocentric latitude of camera measured positive north, degrees 

xp , 6 ,<t> spacecraft azimuth, pitch, and roll Euler angles relative to earth-fixed 

axes, degrees 


Subscripts: 


i 


for camera i = 1 and for object space points i = 2 . . . n + 1 


o 


initial conditions 


S spacecraft 

Superscript: 

T transposing of matrix 

CAMERA ORIENTATION METHOD 

In order to determine the orientation of a camera used to photograph the surface of 
the earth, the following are required: the photographic or image coordinates of at least 
two identification points such as mountain peaks, river bends, and so forth; the earth 
coordinates of these identification points; and the location of the camera with respect to 
the earth. The coordinate system defining both the earth- related identification points 
and the camera is referred to as the object space coordinate system. The photographic 
coordinate system defined by the camera focal axis and the focal plane in which the coor- 
dinates of the images of the identification points may be measured is called the image 
space coordinate system. In order to determine the camera orientation, it is necessary 
to determine the relationship between these two coordinate systems. A detailed discus- 
sion of the mathematical relationships between corresponding points in the image space 
coordinate system and the object space coordinate system is presented in the appendix. 

A least- squares solution to the relationship between the two coordinate systems is also 
shown. The resulting equations can be easily programed into a digital computer system. 

APPLICATION OF METHOD 

The camera orientation method described in the preceding section was applied to 
determine the spacecraft attitude during a portion of the first balloon- launched flight 
test of the Planetary Entry Parachute Program (ref. 4). The test took place at White 
Sands Missile Range, New Mexico. The spacecraft, shown in figure 1, was a 15-ft- 
diameter (4.6 m) 60° half-angle cone. It contained two 16-mm motion-picture cameras 
which are relatively lightweight and inexpensive when compared to other cameras used in 
photogrammetry. Moreover, the cameras were not intended to be used as metric instru- 
ments; instead they were placed on the spacecraft to photograph the inflation character- 
istics of an experimental parachute (ref. 5). Each camera was equipped with a lens 
having a focal length of 10.5 mm, a field of view of 38° by 52°, and a maximum radial 
distortion of less than 2 percent. They were mounted in the base of the spacecraft and 
were pointed toward the ground during the ascending portion of the trajectory. 
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The geodetic coordinates and altitude of the identification points were determined 
by comparing 10- by 7-cm enlargements of the 16-mm frames with a 1:250 000 scale 
contour map. Prominent geographical areas such as mountain ranges and rivers were 
readily identified which established the general region being photographed. Further 
study of the enlargements enabled more minute geographical details such as mountain 
tops and the intersection of rivers to be discerned. 

The relationship between the spacecraft and camera axis systems was fixed prior 
to launch. The cameras were mounted in the base of the aeroshell, pointed aft, and 
rotated 15° from the nominal vertical. 

A total of 170 frames from the two cameras were enlarged and studied for possible 
landmarks which could be used as identification points. After identification of the land- 
marks, the 16-mm film was placed in a film reader which enlarged the film by a factor 
of 20. The coordinates of the identification points with respect to the film (or image 
space) were located and automatically punched on data processing cards. In order to 
increase reading accuracy each point was read three times and the average value of these 
readings was used. 

From two to five landmarks were located and used as identification points on each 
frame investigated. A total of 12 landmarks were identified but only 10 were used in 
this application. The earth coordinates of these points are given in table I. Figure 2 
shows the location of the identification points on a map overlapping the northwest corner 
of the White Sands Missile Range. Enlargements of three of the 16-mm frames analyzed 
with the identification points marked are shown in figure 3. These points are numbered 
and can be compared with the corresponding points of table I and figure 2 . 

The appendix presents the equations which relate the points on earth or object space 
with points on the film or image space. The equations were set up in matrix form and 
solved on a digital computer. The results were then transformed from the camera angle 
orientation system given by o^, a ^d 03 (fig- 4) to the spacecraft attitude system 
given by \p, 9, and cp (fig. 5). 

The relationship between the two orientation systems for this particular application 
is as follows: 

i// = CTj - 180° 

e = -a 2 

(p = -(cr 3 + 15°) 

The data are presented in terms of both the camera orientation system and the spacecraft 
attitude system. 
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The primary purpose of the spacecraft cameras during the PEPP flight test was to 
obtain parachute deployment data and not spacecraft attitude data. Therefore, no attempt 
was made to assure that the focal axes of the cameras were parallel to either the longi- 
tudinal axis of the spacecraft or to each other. Subsequent comparison of the film indi- 
cated that the focal axes of the two cameras were misalined 5° relative to each other. 
After reviewing the assembly procedures and spacecraft drawings, it was estimated that 
the focal axis of each camera could vary as much as 6° with the spacecraft longitudinal 
axis. 


ACCURACY 

A Monte Carlo error analysis was performed to determine the accuracy of the 
camera attitude angles. One hundred frames having two or three observation points were 
forced to be theoretically errorless by adjusting the photographic image coordinates to a 
zero residual condition. (See appendix.) Normally distributed random errors with mean 
zero and expected standard deviation in each of the parameters in the system were intro- 
duced into the theoretical case and new solutions obtained. The sources of error con- 
sidered (and the corresponding deviations) were 

(1) Locating the coordinates of the identification points on the 16-mm film 

(±100 microns) 

(2) Determining the latitude, longitude, and altitude of the identification points 

(±0.5 minutes and ±60 meters) 

(3) Locating the spacecraft by ground-based tracking instrumentation (±60 meters) 

(4) Measuring the focal length of the camera (±1.5 microns) 

(5) Locating the center of the frame (±10 microns) 

Ten cases having random error selection were employed for each frame. The camera 
attitude angle variances were computed by using the following equation: 

V 

Z (°e>ij ' °th,i) 

a„ 2-Jr 1 


where 
i = 1, 2, 3 
77 = 10 
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Histograms of the variances of each angle for the 100 frames were plotted and 
95-percent confidence regions established. The results indicate that aj is known 
within ±1.5°, Og within ±1.0°, and Og within ±3.3°. 

RESULTS AND DISCUSSION 

Data were obtained from both cameras during 28 seconds of the flight test covering 
an altitude range from 130 000 to 143 000 feet (40 to 44 kilometers). Beyond this time 
the spacecraft flight path was nearly horizontal and the cameras were pointed away from 
the surface of the earth. 

Time histories of the camera orientation Euler angles o^, and ag (fig. 4) 
and the aeroshell attitude Euler angles xp, 9, and <p (fig. 5) are presented in figures 6 
to 8. Figure 6 compares the yaw-angle (aj or xp) results obtained from both cameras 
with ground tracking data and a point obtained from an onboard magnetometer. Pitch- 
angle (02 or 0) results are compared with ground tracking data in figure 7 and roll- 
angle (og or (p) data are presented in figure 8. 

At t = 0 the spacecraft was released from a balloon and fell toward earth until a 
cluster of rocket motors in the base of the spacecraft ignited at t = 3.7 seconds. The 
motors burned for about 3 seconds propelling the spacecraft to transonic velocities. At 
t = 7 seconds a parachute was deployed in a reefed condition and remained physically 
connected to the spacecraft until t = 14.2 seconds. After this time the spacecraft was 
in unencumbered free flight and, as shown in figures 6 to 8, its period of oscillation 
increased by a factor of approximately 1.5. 

The Euler angle data from camera 2 agree quite well with the results from cam- 
era 1. Both the amplitude variation during corresponding half- cycles and the oscillation 
frequency are in agreement within the accuracy of this investigation. That the amplitudes 
of the curves at specific times vary by as much as 10° is attributed to the relative mis- 
alinement of the two cameras. A 5° axes misalinement entirely in the Xj - Zj plane 
(fig. 4) of the camera, at a time when the spacecraft pitch angle was 60° and the roll 
angle was near zero, would result in a yaw-angle variation of about 10°. Such a situation 
exists during the time period from 0 to 3.67 seconds as shown in figure 6. 

The camera results are compared with ground tracking data. Radar and optical 
tracking station data were combined to yield flight-path angles y-y and y-g in the yaw 
and pitch ^X E ,Z E ) planes, respectively. These data are represented by the 
dashed lines beginning at 6 seconds in figures 6 and 7. The radar and optical tracking 
stations tracked the spacecraft as a point mass and, therefore, did not obtain the oscil- 
latory attitude variations recorded by the two onboard cameras. Nevertheless, the cam- 
era Euler angles ctj and erg (as well as xp and 6) should oscillate about yy and 
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yp, respectively, in the absence of any wind. A fairing of the mean values of the camera 
data would not coincide with the ground track data. This is attributed to both the offset 
of the camera focal axis from the longitudinal axis of the spacecraft and the lateral veloc- 
ity component of the spacecraft with respect to the ground caused by the wind. The pre- 
vailing winds during the portion of the flight test covered in this investigation were from 
an azimuth between 90° and 110° from north. The spacecraft, according to ground track 
data, flew at a flight-path angle y y of not less than 135°. The resultant velocity vector 
would cause the mean values of and i p to be somewhat north of y y as shown in 
figure 6. In addition, because the winds were horizontal, the mean values of and 6 
would be less than yp during the ascent portion of the trajectory as shown in figure 7. 

CONCLUDING REMARKS 

The use of a relatively unsophisticated onboard camera for postflight determina- 
tion of vehicle attitude within an accuracy of ±3.3° has been demonstrated. The method 
employed requires 

(1) The location of at least two identifiable objects within the view of the camera 

(2) The location of the camera with respect to the objects being viewed 

(3) The orientation of the vehicle with respect to the camera 
The method does not require any other onboard instrumentation. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 18, 1967, 

124-07-03-05-23. 
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APPENDIX 


PHOTOGRAMMETRIC DETERMINATION OF CAMERA ORIENTATION 


In order to adequately describe how the camera orientation can be determined, it 
is necessary to define the coordinate systems in which observations are made. The 
coordinate system in which the position of observed points are known and in which the 
camera is oriented is called the object space. The coordinate system composed of the 
camera focal plane and focal axis and in which image coordinates are measured is the 
image space. (See ref. 6.) The relationship between a point in the object space and 
image space is 


x - x p ** 

y - y P = 


f\ 1 X+ 11 1 Y + i/jZ> 

^X^X + jLlgY + 

/x 2 X + /i 2 Y + 

1X3X + jUgY + 1^3 Z 


(Al) 


The relationship in equations (Al) can be derived with the aid of figure 9. Let Q(X,Y,Z) 
be a point in the object space, then the image of Q will be I (Xj-,Yj,Zjj in the image 
space. Since the origins of the two coordinate systems coincide for all practical pur- 
poses, the following transformation describes the coordinates of Q relative to the 
image space 


X I~ 



x“ 

Y I 

= K 

X 2 I1 2 V 2 

Y 

Z I 


X 3 M3 ^ 3 

Z 


(A2) 


where K is a constant of contraction. In equations (A2) dividing the first and second 
equations by the third equation removes the constant K and gives 


Xj XjX + jUjY + 

Zj ~ X3X + Y + t'gZ 

Yj X 2 X + /i 2 y + ^ 2 Z 
Tj " X 3 X + M3Y + I/gZ 

J 

Letting 

5 " Xp = X I 

y - y P = Y i 

c = Zj 


(A3) 
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and multiplying both sides of equations (A3) by c gives 


x - Xp = c 


y - y P 


= C 


A x x+ Y+ V) Z\ 

yXgX + /igY + I/gZy 

(\^X. + m 2 Y + 
i XgX + MgY + t*gZ , 


(A4) 


Equations (A4) are identical to equations (Al). 

For each known point in the object space there exist two equations (Al) relating 
the object space to the image space. The directional cosines Xj, n^, v j, . . ., of 
the image space axes relative to the object space axes can be expressed in terms of 
three angles a y Og, and Og. Consequently, these angles can be used to describe the 
orientation of the camera coordinate system relative to the object space coordinate sys- 
tem. To obtain the directional cosines consider three successive rotations through the 
angles o^, o^, anc * °g- Using figure 4 and imposing a constant of contraction for the 

camera, the following transformation is obtained: 


- — 1 
X - Xp 


-COjCOg - SOjSOVjSCTg 

SOjCOg - COjSOgSCTg 

Co 2 Sa 3 


r n 

X 


^l^l^l 

X 

y - y P 

= K 

Ca^SOg - So^SCTgCo-g 

-So^SOg - Ccr^SOgCag 

Ccr 2 Co 3 


Y 

= K 

x 2^2 v 2 

Y 

c 


SOjCOg 

COjCOg 

S°2 


Z 


XgMg^g 

Z 


where, for convenience, C and S represent cosine and sine, respectively. This equa- 
tion is equivalent to equation (A2) . 


The Least- Squares Solution 

The solution of equations (Al) contains six parameters which under theoretical con- 
ditions are constants. These parameters are Xp, yp, c, o^, Og, and Og of which 
Xp, yp, and c are measured independently for this experiment and are not unknowns 
in the solution. A solution of equations (Al) for <7j, Og, and Og can be found with two 
properly chosen observations. Associated with x and y in equations (Al) are errors 
6 and e. Since these errors exist, a computational method is needed which yields the 
best possible results with all the information available. The method of least squares 
which is described subsequently uses a minimum error criterion and has been used in 
the data reduction for this investigation. 
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In general, equations (Al) with the associated errors can be written as 

(A 5) 

where F and F are nonlinear functions of Oj and in order to find a solution they 
must be linearized. Expanding equations (A5) in a Taylor's series about a nominal set 
Uj Q and dropping the higher order terms results in the following linear approximations: 


x i = F ( a j) + e i 
yi = F ( CT i) + | i 



Ax i = x i - x i,o = b il( CT l - \o) + b i2 ( °2 - a 2,o) + b i3 ( °3 " °3,o) + e i ? 
A ^i = h ~ h,o = b il( CT l - a l,o) + b i2 ( ct 2 - ct 2,o) + b i3( CT 3 " a 3,o) + «i| 


(A6) 


where 


b 


9F, 


il SCTj 


a, cr 0 CT„ 
l,o 2,0 3,o 


b. * = 


8F 4 


il 9m 


u i o o 
l,o 2,o o,o 


dF i 

b.„ = — - 
12 9CTo 


<J-i CT„ 0T o 

1,0 2,0 3,o 


9F. 


b. 0 = 


i2 d(7 n 


a l,o a 2,o%o 


b i3 = 


9Fi 


9ao 


a a 0 a„ 
l,o 2,o 3,o 


9Fi 
b i3 = eon 


a- a 0 a Q 
l,o 2,o 3,o 


Letting 

A<T i * (°i - a i,o) 

equations (A6) can be put in the following form: 


Axi = ^ byAOj + Ci 

3=1 

3 

Ayi = I b ij Aa j + ? i 
3=1 


(i = 1, . . ., n) 


> (A7) 


For further considerations the linear equations (A7) corresponding to the ith observation 
are expressed in matrix notation where 
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Vj = BjAcx + (A8) 



Axj 


b il 

b i2 

b i3 

Vi = 


B i = 




AYi 

!_ _l 

|_ b il b i2 

b i3 


"ACTj 








fe i' 



Act = 

Act 2 

e i = 

Ji_ 




Act 3 





Then for n observations there are n matrix equations of the form of equation (A8) 
which may be written 

V = BAct + e (A9) 

where 


' v i n 


Bl 


e i 

v 2 


b 2 


e 2 

• 

B = 

* 

e = 

• 

v n 

i 

B n 


_ e n_ 


The problem may be restated: given V and B find the best estimate A o for Ac, 

The best estimate Act is the value of Act which minimizes the sum of the squares 

- T- 

of the residuals e e where 

e T e = (v - BAct) T (v - BAct) (A10) 


In order to minimize equation (A10) the first variation 6 with respect to Act must 
vanish, that is 


6(e T e) = s[(v - BAct) T (v - BAct)_ 
6(e T e) = -2 (v T - Act T B T )b5Act = 0 


= 0 


(All) 


Since 5 Act * 0, equations (All) can be satisfied if 
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or 

B T BAct = B T V (A12) 

Solving for the estimate of A a in equation (A12) gives 

Act = (b T b) 1 B T V (A13) 

A second necessary condition for equation (A10) to be a minimum is that the second 
variation with respect to A a be positive definite. Upon examination the second varia- 
tion is 

6 2 (e T e) = 26 Act T B T B5Act 

which is positive definite. Therefore, equation (A13) is a valid expression for Aa. 

Since equation (A13) is based on a linear approximation with nominal o^ 0 , Aa 
can be used to find the best estimates by With the relationship a = a 0 + A a, the value 
of Aa which minimized equation (A10) leads to a new nominal a^ Q = a-^ 0 + Aa. This 
process implies an iterative procedure which continues until Aa 0 and the value of 
aj j0 that leads to this result is the best estimate of &j for aj. 


Partial Derivatives of Projection Equations 
For equations (A5) 


*i = F (°j) + e i 

h = F (°j) + e i 


(j = 1, 2, 3) 
(i = 1, . . n) 


The partial derivatives of F and F with respect to a-^, a 2 , and ag are as follows. 

Let 

p = X^X + jjl^Y + v^L 
q = X 2 X + fi 2 Y + 
r = XgX + jUg Y + Z/gZ 


then 


9F 

da. 


a P 

9ai 


r - 


3r 

3ai 


P 


1 r 2 

/ 3p 9r 

c ( — =— r p 

8F _ \ aa 2 da 2 

9a 0 t»2 
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8F 

8o 3 



8r 

8o 3 


9F 

3a^ 



8F 

8ct 2 



8F 

9c 3 



w x = "i x - x i Y 

3q 

^-=m 2 x-x 2 y 


■|3L = MqX - XoY 
8o^ 3 3 




8ao 


z^X sin cfj + Y cos ct-jJ - (sin o 2 sin 


•|SL = _^ 2 (x sin CTj + Y cos OjJ - (sin cr 2 cos 


= -^(X sin + Y cos a-j^j + Z cos a 2 

2 

lij = X 2 X + "2 Y + "2 Z = N 

^=-(X lX+)ll T + ^Z).-M 
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Coordinate Transformation 

Object space points and the position of the camera are initially identified in terms 
of geodetic latitude, longitude, and altitude above sea level. These data are obtained from 
maps of the photographed area and radar observations of the vehicle trajectory. In 
order to reference the data relative to the camera as described previously the earth 
centered geocentric coordinates of both the object space points and the vehicle position 
are computed by 


Xi = (Ri + h^cos cos Aj 
Yi = (Ri + h^cos <f>! sin A* 

Zi = (r^ + hj^sin 

and 

- ll’35.6635" sin 2^ + 1.1731” sin 4$i - 0.0025" sin 6<f>i 
Ri = 6378.388(0.998320047 + 0.001683494 cos 2$i - 0.000003549 cos 4*) 


where 

i = 1 for the vehicle position 

i = 2 . . . n+1 for the object space points 

In the preceding the X axis is in the equatorial plane pointing toward the 
Greenwich meridian, the Y axis 90° east in the equatorial plane and the Z axis 
toward the north pole. A slight error is introduced by adding altitude to the earth's 
radius vector but this is negligible for the accuracy desired of this system. 

The following transformation maps the geocentric coordinates into object space 
coordinates relative to the camera where the X axis points east, the Y axis north, 
and the Z axis toward the zenith. 


x,r 


-sin Aj cos Ai 0 


x + i - *r 

Y i + i 


-cos Aj sin -sin Aj sin cos 


Y i + i - y i 


! 

cos Aj cos <f>i sin Aj cos sin 


3 + i - \ 


This transformation is obtained by a positive rotation about the Z axis through 
an angle of 90° + Aj followed by a positive rotation about the new X axis through an 
angle of 90° - 4>j. 
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TABLE I.- GEODETIC COORDINATES AND ALTITUDE 


OF IDENTIFICATION POINTS 
[See fig. 2j 


Point 

Latitude 

Longitude 

Altitude 

ft 

m 

2 

33°13.8' 

106°26.5' 

4050 

1234 

4 

33°17.9' 

106°31.9' 

8358 

2548 

5 

33°16.8' 

106°23.6' 

4075 

1242 

6 

33°40.3' 

106°58.7' 

4912 

1497 

7 

33°29.0' 

107° 7.8’ 

4500 

1372 

8 

33° 9.0’ 

107°11.0' 

4500 

1372 

9 

33°13.8' 

106°38.4' 

7600 

2316 

10 

33°20.0’ 

106°30.5' 

7900 

2408 

11 

33°36.6' 

106°21.3' 

8114 

2473 

12 

34°15.3' 

106°53.9' 

4750 

1448 
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